library(ez)
options(contrasts = c("contr.sum","contr.poly"))
score<-c(1,2,3,4,5,6,8,7,9,10,12,11,13,14,15,16,17,18,20,19,21,22,24,23)
A <- c(rep(1,12),rep(2,12))
B <- rep(c(1,1,1,1,2,2,2,2,3,3,3,3),2)
subject <- c(rep(c(1,2,3,4),3),rep(c(5,6,7,8),3))
rm.data<-data.frame(subject,A,B,score)
rm.data
rm.data$A<-factor(rm.data$A)
rm.data$B<-factor(rm.data$B)
rm.data$subject<-factor(rm.data$subject)
out2 <- ezANOVA(data=rm.data,dv=.(score),wid=.(subject),within=.(B),between=.(A),type="II")
out3 <- ezANOVA(data=rm.data,dv=.(score),wid=.(subject),within=.(B),between=.(A),type="III")

Table1505 <- read.csv("http://www.statpower.net/Content/311/Lecture Notes/Table1505.csv")
rm.data <- data.frame(melt(Table1505,id.vars=1:2,measure.vars=3:6))
rm.data$A <- factor(rm.data$A)
colnames(rm.data) = c("Subject","Method","Time","Score")
ezANOVA(rm.data,wid=.(Subject),dv=.(Score),between=.(Method),within=.(Time))